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RAYLEIGH-BRILLOUIN  SCATTERING  AND  STRUCTURAL 
RELAXATION  OF  A VISCOELASTIC  LIQUID. 


Y.-H.  Lin  and  C.  H.  Wang 


Department  of  Chemistry 
University  of  Utah 
Salt  Lake  City,  Utah  84112 


Abstract 


The  linear  response  theory  using  a complete  set  of  dynamic  variables 
involving  density,  velocity  and  energy  fluctuations  is  used  to  analyze  the 
Rayleigh-Brillouin  spectrum  of  a viscoelastic  liquid.  An  exact  equation  in- 
cluding the  effect  of  structural  relaxation  has  been  obtained.  The  result 
is  discussed  in  the  fast  and  slow  relaxation  limits.  In  the  former  case, 
the  Rayleigh-Brillouin  spectrum  is  identical  to  that  predicted  by  the  classi- 
cal hydrodynamic  equations,  whereas  in  the  latter  case  a new  structural  cen- 
tral peak  is  found,  in  addition  to  the  anomalous  dispersion  and  relaxation 
effects  present  in  the  frequency  and  linewidth  data  associated  with  the  Bril- 
louin  peak.  The  evolution  of  structural  relaxation  and  its  effect  on  the  en- 
tire Rayleigh-Brillouin  spectrum  is  described.  The  structural  central  peak 
is  most  pronounced  with  the  frequency  dispersion  and  the  linewidth  maximum 
are  present.  The  theoretical  result  has  been  used  to  calculate  the  Rayleigh- 
Brillouin  spectra  of  polypropylene  glycol  at  various  temperatures.  The  re- 
sults of  the  frequency  shift,  the  spectral  linewidth  and  the  Landau-Placzek 
ratio  are  in  good  agreement  with  the  experiment. 


' 7 TFS  j 

cSTl 


I.  Introduction 


The  study  of  molecular  dynamics  of  a viscoelastic  liquid  using 
the  Rayleigh-Brillouin  scattering  technique  is  of  current  research 
interest.  Rayleigh-Brillouin  scattering  arises  from  density  fluctu- 
ations and  the  scattering  spectrum  contains  information  about  the 
translational  motion  of  an  ensemble  of  molecules.  In  a viscous  molec- 
ular liquid  consisting  of  polymeric  molecules,  interactions  between 
molecules  as  well  as  between  segments  of  the  same  molecule  play  a very 
important  role  in  the  Rayleigh-Brillouin  spectrum.  In  such  a system, 
a considerable  frequency  dispersion  as  well  as  a maximum  in  the  line- 
width  of  the  Brillouin  doublet  are  known  to  occur.  Accompanied  with 
the  frequency  dispersion  and  linewidth  anomalv,  one  also  finds  sig- 
nificantly more  intensity  distributed  in  the  region  between  the  un- 
shifted Rayleigh  and  Brillouin  peaks  than  that  found  in  a simple 

molecular  liquid.  These  experimental  results  cannot  be  accounted  for 

2 

in  terms  of  the  classical  hydrodynamic  theory  of  simple  fluids. 

To  explain  the  anomalous  sound  dispersion  and  attenuation  maximum 

3 

of  a highly  viscous  liquid,  Isakovich  and  Chaban  have  provided  a phe- 
nomenological theory  in  which  the  viscous  liquid  is  modeled  as  a "micro- 
inhomogeneous  medium"  consisting  of  several  regions.  The  occurrence  of 
diffusion  and  energy  exchange  mechanisms  between  different  regions  of  the 
medium  is  shown  to  be  responsible  for  the  attenuation  of  the  hypersonic 

wave.  Although  this  theory  can  qualitatively  account  for  the  dispersion 

4 

and  attenuation  phenomena  of  the  hypersonic  wave  in  polymeric  liquids, 
it  does  not  have  a satisfactory  statistical  mechanics  foundation. 

We  have  recently  used  the  linear  response  theory  to  analyze  the 
Brillouin  scattering  results  of  bulk  liquid  polypropylene  glycols  (PPG) 


2 


and  polyethylene  glycols  (PEG).  The  analysis  was  carried  out  using  a 
two  coupled  variables  (density  and  velocity)  theory,  and  the  theoreti- 
cal result  was  found  capable  of  giving  a satisfactory  interpretation 
of  the  frequency  dispersion  and  linewidth  data  in  both  PPG^  and  PEG.*5 
It  was  shown  that  the  second  moment  of  the  intersegraental  interaction 
potential  and  its  associated  relaxation  dynamics  are  responsible  for 
the  anomaly  found  in  the  linewidth  and  frequency  shift  of  the  Brillouin 
peak.  The  theoretical  results  cannot,  however,  account  for  the  entire 
Rayleigh-Brillouin  spectral  distribution,  due  to  the  fact  that  energy 
fluctuation  is  not  included  in  the  theoretical  analysis. 

In  this  paper,  we  report  the  result  of  a complete  theoretical  anal- 
ysis of  the  Rayleigh-Brillouin  scattering  spectrum  using  the  linear  re- 
sponse theory.  This  theoretical  result  includes  the  energy  fluctuation 
and  is  applicable  not  only  to  the  frequency  shift  and  linewidth  data  of 
a polymeric  liquid  but  also  to  the  intensity  ratio  between  the  Rayleigh 
and  Brillouin  components.  We  also  show  how  the  structural  relaxation 
mechanism  can  give  rise  to  a new  central  peak.  This  new  central  peak 
is  due  to  coupling  between  the  translation  and  the  structural  rearrange- 
ment motion  in  a viscoelastic  liquid.  We  draw  an  analogy  between  this 
dynamic  structural  central  peak  and  the  central  peak  associated  with  the 
structural  phase  transition  in  ferroelectrics  and  the  translation  rota- 
tion coupling  in  orlentatlonally  disordered  solids.  We  show  how  the 
central  peak  disappears  in  the  limit  of  fast  modulation. 

II . Theory 

For  a viscoelastic  liquid  made  up  of  a collection  of  identical 
scattering  elements,  the  isotropic  part  of  the  Rayleigh-Brillouin  spec- 
trum observed  at  scattering  angle  9 is  determined  by  the  qth  mode  of 


r 

* 


the  Fourier  transform  of  the  time  correlation  function  of  density  fluc- 
tuations given  by 

C (q , t)  - <6p(q,t)5p  (q,o)>  (l) 

PO  - 

where  6p(q,t)  is  the  qth  mode  of  the  fluctuation  of  density,  q being 
the  scattering  vector  whose  amplitude  is  given  by  |q|  =»  sin  -r-. 

Here  is  the  wavelength  of  the  incident  light  and  n being  the  index 
of  refraction  of  the  liquid.  In  this  section,  we  shall  studv  the 
density  fluctuation  in  viscoelastic  liquids  by  calculating  the  qth 
mode  of  the  densitv  fluctuation  spectrum. 

For  an  isotropic  fluid  svstem,  it  is  easy  to  see  bv  svmmetrv  con- 
siderations that  density  (5p),  velocity  (V)  and  energy  (5c)  fluctuations 
form  a set  of  three  dynamic  variables,  which  do  not  couple  to  any  other 
variables.  These  are  given  respectively  by 

5p  - T.  eiqzj  (2) 


V - Z z ,eiqZ j 

j J 

6e  - Z Psm ,(x  2 + y.;  + z 2)  + H E U11>eiqZj 
a J J J J i+a 


where  in  Eqs.  (2)-(4)  the  direction  of  the  scattering  vector  q is  chosen 
to  be  along  the  z-direction.  , y^  and  z^  are  the  coordinates  of  the 
jth  scatterer.  In  a polymeric  fluid  these  can  be  regarded  as  the  coor- 
dinates of  the  jth  segment  of  a polymer  chain.  V is  the  component  of 
velocity  in  the  direction  of  q,  as  the  transverse  components  will  not 
be  coupled  to  the  density  fluctuation  and  thus  have  no  effect  at  all  on 
the  polarized  Rayleigh-Brillouin  spectrum  of  an  isotropic  fluid  system.7 

is  the  interaction  potential  between  scatterers  i and  j.  In  Eqs. 
(2)-(4),  the  q dependence  in  6p,  V and  6c  is  not  written  out.  explicitly. 


The  sec  of  variables  defined  in  Eqs.  (2)-(4)  (*5p,  V,  5e)  gives  rise 

Co  che  collecCive  mocion  of  a sound  wave.  The  dynamic  behavior  of  Che 

collecCive  variables  are  mosdy  conveniencly  described  by  relaxacion 

funcCions  decermined  by  variables  noc  included  in  <5p,  V,  and  5e.  To 

analyze  Che  effecc  of  relaxacion  on  6p,  V and  Se,  we  use  Che  Zwanzig- 
0 

Mori  formalism.  We  shall  show  chac  Che  final  equacions  obcained  are 
capable  of  describing  Che  sound  dispersion  and  accenuacion  behavior 
in  a polymeric  liquid.  We  shall  also  show  chac  in  Che  limic  of  fast 
relaxation  the  equations  reduce  Co  Che  usual  hydrodynamic  resulcs. 


A)  Dynamic  equacions 

We  start  from  the  assumption  that  for  the  time  scale  of  interest, 
Che  dynamics  of  the  system  is  governed  by  the  evolution  of  the  set  of 
secular  variables  {6p,  V,  $e}.  We  use  the  Zwanzig-Mori  projection 
operator  cechnique  Co  calculate  the  motion  of  secular  variables.  Be- 
fore proceeding,  we  note  that  the  quantity  <6eSp**  is  non-zero.  There- 
fore, we  use,  instead  of  6p,  V and  5e,  a set  of  orthonormal  variables 


A given  by 


where 


- 5p/  <| 5p  I2*** 

£2  - v/  <|  v|  2>4 


£3  - 77 -q  - n^] 

< I vSc  I * 


and 


<5 z6o*> 


< ! 3 e | ->"*  < ! 5p  I 


(9) 


From  the  definition  given  above,  one  has  <A  A > =*  I,  where  I is 
a 3x3  unit  matrix.  The  angular  brackets  indicate  an  ensemble  average. 
We  show  in  Appendix  I that  if  local  equilibrium  is  assumed  £3  is  pro- 
portional to  the  temperature  fluctuation. 

For  the  column  vector  A of  dynamic  variables,  the  equation  of 
motion  is  given  formally  by  the  Liouville  equation 


A(t)  - iLA(t) 


(10) 


where  L is  the  classical  Liouville  operator  governing  the  time 
evolution  of  A(t).  While  it  is  straightforward  to  write  down  a 
formal  expression  for  L in  terms  of  the  kinetic  and  potential  energy 
operators  of  the  system,  in  practice  this  is  seldom  done  for  a com- 
plex molecular  system  such  a3  a polymer,  due  to  its  large  number  of 
internal  degrees  of  freedom  and  the  interraolecular  and  intramolecular 
potential  functions.  For  the  present  work  it  is  only  necessary  to 
utilize  the  formalistic  symmetry  properties  of  the  Liouville  operator 
such  as  time  reversal,  reflection  and  translation  symmetry.  In  its 
formal  definition  the  Liouville  operator  is  self-adjoint. 

Following  the  well  known  procedure  the  equation  of  motion  for 

8 

A(t)  can  be  shown  to  follow  a generalized  Langevin  equation. 


~ A(t)  - i Q A(t) 


dr  K(t)*A(t-r)  + F(t) 


(ID 


1 


where  the  frequency  matrix  is  given  by 


6 


Q = < (LA) A > , 

other  quantities  are  given  by 


(12) 


K(x)  - <F(t)F+> 

(13) 

F(t)  > eiQLTF 

(14) 

F - iQLA 

(15) 

where  Q - (1-P),  P being  the  projection  operator  defined  as 


PG  3 <GA  > <AA  > A 


In  the  small  q limit,  Eq.  (14)  can  be  approximated  by 


F(t)  - eiLT  F 


(16) 


(17) 


Substituting  Eq.  (5)  into  Eq.  (12),  and  with  the  help  of  Eqs.  (6)-(8) , 
we  obtain  the  frequency  matrix 


/ 0 UT  0 \ 


“T 


0 s 


S*  0 


(18) 


where 


W — w 

/ NkT | , , 2 ^ 

“t " ) “I50!  * 


(19) 


Equation  (19)  reduces  to  the  well  knovm  result  q(pm\^,)  in  the 
small  q limit.  Here  Is  isothermal  compressibility. 

The  parameter  S is  given  by 
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S - <(U2K3*> 


<(LV)<5e*> 

Z U "nu)T 

< | V | 2>^<  j oe | T 


and  is  real  only  in  small  q limit.  We  shall  show  shortly  that  S 
is  related  to  the  difference  between  the  isothermal  and  adiabatic 
sound  frequencies. 

The  dissipation  part  can  be  calculated  by  the  following.  Accord- 
ing to  Eq.  (15),  one  can  readily  find  that 


where 


F - QiLA  =*  1LA  - iOA 


f2  ” iL£2  “ ~ is?3 


f3  = iU3  - iS*£, 


Thus  the  memory  function  matrix  K(t)  has  the  general  form. 


k22(t)  0 
0 K33(t) 


where  K22(t)  - <f2(r)f2*>  and  K33(t)  - <f3(t)f3*>.  K23(t)  and  K32(t) 

vanish  at  t*0  and  are  assumed  to  be  negligible  at  finite  r. 

Equation  (11)  together  with  the  frequency  matrix  given  in  Eq. 
(18)  and  the  memory  function  matrix  given  in  Eq.  (24)  can  be  solved 
by  the  Laplace  transform  technique.  Taking  the  Laplace  transform 
on  both  sides  of  Eq.  (11),  we  obtain, 


X(z)A(z)  - [zl  - in  + £(z)l  A(z)  - A(o)  + F(z) 


where  A(z)  and  K(z)  are  the  Laplace  transforms  of  A(c)  and  K(t), 
respectively.  The  expression  for  X(z)  is  given  by 
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■ ioj^,  z+K22  (z)  -iS 


a+K33 (z) 


The  Rayleigh-Brillouin  spectral  density  is  equal  to  the  real  part 
of  '•Ci  ( -^ ) C 1 * > • Thus  we  only  need  calculate  this  quantity  from 

Eq.  (25).  From  Eq.  (25)  and  also  using  the  fact  that  <F(zHi*>  * 0, 
we  obtain  after  some  calculation. 


<Cl  (z)Ci*>  = "t  v r {(z+K22(z))(z+K.33(z))  + S2} 


where  ; X I is  the  determinant  of  X(z),  and  is  given  by 


Xj  * [z(z  +K?0(z))  +u  2]  [z  +K33(z)]  + zS‘ 


The  Rayleigh-Brillouin  spectral  density  is  thus  equal  to 


I(q,ui)  * Re<?i  (z)£i*> 


[z  + K22(z)l  [z  +K33(z)1  + S2 

- Re  (S 

Iz(z+K22(z))  + <»T2]  (z  +K33(z)1  + zS2!.,^ 

Equation  (29)  is  the  dynamic  equation  describing  the  density- 
density  fluctuation.  This  equation  will  be  used  to  interpret  the 
Rayleigh-Brillouin  spectra  of  a viscoelastic  liquid  in  Section  IV. 


B)  Adiabatic  frequency 

Before  proceeding  to  study  the  property  of  the  dynamic  equation, 
we  want  to  associate  the  quantity  S with  the  well  known  physical 
quantity.  One  can  relate  the  quantity  S to  adiabatic  sound  frequency 


I 


by  diagonalizing  Eq.  (18).  This  gives  rise  to  three  eigenvalues 


X+  =»  + S~  « u>s 


X * — + S2  * -u)g 


Xn  - 0 


(30) 

(31) 

(32) 


The  corresponding  normalized  eigenvectors  are  found  to  be 


(33) 


and 


/ - \ 

■ a)  1 


S_ 

Js 

0 


\ 


(34) 


T / 


In  the  absence  of  relaxation  effects,  Sj  can  be  expressed  in 


terms  of  and  $ as 
± o 


S,  ^ — r + ) 

<|6p|2>%  /T  ' “ S + “s  ws  °l 


(35) 


The  eigenvectors  <f>+  and  $ are  associated  with  undamped  acoustic 
waves  propagated  with  frequency  wg  and  -uig  respectively.  The  third 
eigenvector  is  associated  with  the  unshifted  Rayleigh  central  peak. 
The  intensity  of  each  peak  is  proportional  to  the  square  of  the  co- 
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efficient  of  the  corresponding  eigenvector.  Thus,  with  the  help  of 
Eq.  (35),  we  obtain  the  Landau-Placzek  ratio  as 


LP 


1 

c_ 

2 Ir 


S2 


V 


- i 


(36) 


Since  the  Landau-Placzek  ratio  in  the  absence  of  relaxation  is 

Q 

equal  to  y-1,  where  y is  p/Cy,  it  follows  immediately  the  result 
ojs  = /ywT,  which  is  the  adiabatic  sound  frequency. 


III.  Effect  of  Relaxation 

To  consider  the  effect  of  relaxation  we  need  return  to  Eq.  (29). 
The  resonance  property  of  I(q,ai)  depends  significantly  on  the  quanti- 
ties K22 (itu)  and  K3  3(iw)  as  well  as  the  values  of  14^  and  u)  . This 
can  be  seen  by  equating  to  zero  the  denominator  on  the  right  hand  side 
of  Eq.  (29)  after  setting  z ■ iw. 


(iu>)3  + (iw)2  1X2(10)  +X3C0J)]  + iu>  [ X 2 (to)  X 3 (w)  + u>s2]  + 


to^j,2X 3 (to)  * 0 


(37) 


where 

Xi(w)  " Kii(:U)  *["  e_iU>t  O8) 

Because  of  Eq.  (38),  Eq.  (37)  is  in  general  a complex  equation 
of  frequency  and  does  not  yield  a simple  solution.  Nevertheless,  it 
can  be  shown  that  in  the  limit  of  small  q,  the  memory  function  is  a 
decaying  function  of  time  and  for  simplicity  can  be  approximated  as 

Kii(t)  " v'*^1 

where  - K^o),  and  are  the  correlation  times  given  by 


(39) 
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:i-iC77  r Kii<t>dt 

ii(o)  'o 


In  this  approximation,  we  obtain  from  Eq.  (38) 


i i 1 +iuixi  ' ' 

Using  Eq.  (41),  one  sees  that  Eq.  (37)  is  at  least  fourth  order 
in  and  it  leads  to  a four  peak  structure  in  the  Rayleigh-Brillouin 
spectrum.  We  now  consider  two  limiting  cases. 


A)  Fast  relaxation,  Vx^>w-< 


In  this  case  becomes  a^x^  and  Eq.  (37)  reduces  to 


(iw)3  + (iu>)2  [02X2  + 013X3]  + (ito)  [<*2X2013X3  + wg2l  + 


“r  <*3T  3 - 0 


which  is  exactly  the  results  obtained  from  the  linearized  equations  of 
2 

hydrodynamics.  Rayleigh-Brillouin  scattering  spectra  predicted  by 

« , 7 

this  equation  have  been  discussed  extensively  in  the  literature. 

Equation  (42)  has  three  roots:  two  correspond  to  the  Brillouin 

side  bands  and  one  to  the  unshifted  Rayleigh  peak.  For  most  liquids  the 

condition  u>  » o,x.  (i  • 2,3)  applies,  and  one  can  obtain  the  three 
s i i 

roots  of  Eq.  (42)  by  a perturbation  technique.  The  results  are  given 


“ ■ ±“  + i j C C 1 — ~ ) 03x3  + a2X2l 


a)  » + i—  03X3 

o y 


Thus  we  note  that  in  the  fast  relaxation  limit,  due  to  the  large 
value  of  tyx^  (i  - 2 or  3)  the  peak  frequency  of  the  sound  wave  is  de- 


a 
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termined  by  >ag  and  its  width  by  a^T^.  The  width  °i  the  unshifted 
Rayleigh  peak  is  determined  only  by  ot 3T  3 / 

T*  7 

To  give  a comparison  with  the  hydrodynamic  result,  one  makes 
the  following  identifications: 


where  Dv  and  D^,  are  the  longitudinal  kinematic  viscosity  and  the 
thermal  diffusivity,  respectively.  These  are  related  to  shear 
viscosity  ( ng),  bulk  viscosity  ( n ) and  thermal  conductivity  (A)  by 


" (nv  + 3 ns)/mpo 


D_  - A/m  p c 
l op 


(48) 


Here  the  various  symbols  have  their  usual  meanings. 

Using  the  expressions  for  K^(t)  (Eqs.  (22)  and  (23))  one  can 

also  obtain  microscopic  expressions  for  the  transport  coefficients. 

The  microscopic  expressions  for  the  transport  coefficients  have  been 

given  in  the  literature.  The  dynamics  of  the  density  fluctuation  in 

9 

the  fast  relaxation  limit  was  first  discussed  by  Mori. 


B)  Slow  relaxation,  Vt2<uj-w3 

This  is  the  relevant  case  for  most  viscous  liquids,  and  under- 
standing of  this  case  is  necessary  for  analysis  of  the  Rayleigh- 
Brlllouin  data  in  polymeric  liquids. 

Since  the  dynamics  of  thermal  conductivity  processes  has  only  a 
negligible  effect  on  the  relaxation  behavior  of  the  sound  wave,  the 
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frequency  dependence  of  A3  can  be  ignored.  Thus  one  can  set 


1J 


l3  * a 3T 3 

a2T2 
k 2 1+iw  T , 


where  x,  will  be  referred  to  as  the  structural  relaxation  time. 

One  further  substitutes  Eqs.  (49)  and  (50)  into  Eq.  (37)  to  find  its 
solution  for  w.  The  computer  calculation  of  the  spectrum  using  the 
exact  expression,  Eq.  (29),  has  shown  that  (i=2,3)  is  only  a 
fraction  of  u>g2  for  a polymeric  liquid  (see  Section  IV).  Therefore, 

Eq.  (37)  can  be  solved  as  a perturbation  of  ct^/u>g2.  The  perturbation 
solution  of  Eq.  (37)  to  first  order  in  ci^/u^2  ( i =2,3)  gives  four  roots 

u>l  * ^ a3T 3 “ iDT92  (51) 

+ 0t2a)sT22  i 

U * ~ 2 [ 1 + (u)gx2)2]  + 2 [ 1 + (wsx2)2] 


{a2x2  +(!--)  (a3x3  +us2x22ci3X3)} 


“2  * ^2  " l+wg2x22 


One  sees  that  the  root  aij  is  identical  to  the  unshifted  Rayleigh 
peak  given  in  Eq.  (44),  indicating  that  it  is  free  from  the  effect  of 
structural  relaxation  provided  that  the  frequency  dependence  of  X3  is 
ignored.  The  roots  correspond  to  a high  frequency  doublet  with 
resonance  centered  around  the  real  part  of  10^  and  the  width  given  by 
the  imaginary  part.  The  peak  frequency  as  well  as  the  linewidth  of 
the  Brillouin  doublet  are  modified  significantly  by  the  effect  of 


E 
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structural  relaxation.  One  notes  that  Eq.  (52)  reduces  to  Eq.  (43) 

in  the  limit  of  ui  t0-»o. 

s 2 

A new  feature  in  this  case  is  the  presence  of  an  additional 

central  peak  with  its  width  given  by  Eq.  (53).  The  width  of  this 

central  peak  is  determined  by  the  structural  relaxation  and  shows 

a frequency  dependence  due  to  its  coupling  to  the  sound  wave. 

The  perturbation  solution  given  in  Eqs.  (51)— (53)  are  new  results. 

These  are  important  for  the  interpretation  of  the  Brillouin  scattering 

spectra  of  viscoelastic  liquids.  We  should  remark  that  these  roots 

differ  from  Mountain's  results  on  thermal  relaxation  arising  from  the 

coupling  of  internal  degrees  of  freedom  of  molecules  to  transla- 
10-12 

tion.  *■  Mountain  has  shown  that  such  a coupling  gives  rise  to 
a relaxation  central  peak  whose  linewidth  is  given  by  (v^/v)2l/r, 
where  and  v are  the  low  frequency  adiabatic  sound  velocity  and 
the  sound  velocity  at  finite  frequency,  respectively;  t is  the  thermal 
relaxation  time.  Our  expression  is  somewhat  simpler  than  Mountain's 
expression  for  thermal  relaxation. 

The  Brillouin  frequency  and  linewidth  given  in  Eq.  (52)  also 
differ  from  Mountain's  result.  Despite  the  fact  that  Mountain  deals 
with  thermal  relaxation  due  to  coupling  of  the  internal  degrees  of 
freedom  to  translation  and  we  are  concerned  with  structural  relaxation 
due  to  coupling  of  internal  as  well  as  external  degrees  of  freedom 
to  translation,  the  difference  between  the  two  results  is  due  to 
different  methods  used.  In  the  present  work  a general  statistical 
mechanical  theory  is  employed  to  analyze  the  structural  relaxation 
phenomena,  whereas  in  Mountain's  calculation  a modified  Navier-S tokes 
equation  allowing  for  a frequency  dependent  bulk  viscosity  is  used. 


I 


" I : 


The  validity  of  the  modified  hydrodynamic  equation  in  a form  given 
in  ref.  10  is  yet  to  be  established. 
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It  is  interesting  to  note  that  structural  relaxation  has  a 
significant  effect  on  the  Rayleigh-Brillouin  spectrum  of  a polymer 
liquid  consisting  of  flexible  polymer  chains.  At  high  temperature 
the  structural  relaxation  rate  is  fast  compared  with  the  frequency  of 
the  hypersonic  wave  and  in  this  case  the  structural  relaxation  con- 
tributes mainly  to  the  Brillouin  doublet  (Eq.  (5?))  through  the 
02T2  term.  In  this  case  the  intensity  of  the  structural  central 
peak  is  smeared  out  over  a wide  range  of  frequency  and  is  not  detect- 
able (see  Fig.  3A) . At  intermediate  temperature  when  the  structural 
relaxation  rate  l/x„  is  the  order  of  u the  linewidth  of  the  Brillouin 
component  doublet  displays  a maximum.  This  can  be  demonstrated  by 
taking  the  derivative  of  the  imaginary  part  of  Eq.  (52)  with  respect 
to  x 2 and  then  setting  the  result  to  zero.  This  gives  in  terms 
of  o . In  this  case,  the  hypersonic  velocity  displays  considerable 
dispersion  (c.f.  Eq.  (52)],  and  the  linewidth  of  the  structural 
central  peak  is  approximately  equal  to  Wg-  Thus  the  structural  central 
mode  contributes  significantly  to  the  intensity  of  the  Brillouin 
doublet  (Fig.  3B)  as  well  as  the  intensity  in  the  region  between  the 
unshifted  Rayleigh  and  the  Brillouin  peaks.  At  low  temperature,  on 
the  other  hand,  the  structural  relaxation  rate  is  slow  compared  to  the 
hypersonic  frequency,  but  fast  compared  to  the  thermal  diffusion  rate; 
in  this  case  the  structural  relaxation  manifests  itself  as  a broad  line 


centered  at  the  laser  frequency  (see  Fig.  3D).  Lastly  at  temperatures 
close  to  the  glass  transition  temperature,  the  structural  relaxation 
is  very  slow  or  nearly  frozen-in,  and  the  structural  central  peak 
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becomes  very  narrow  and  hidden  under  the  ordinary  Rayleigh  diffusion 
peak  (see  Fig.  3F) . 

We  remark  finally  that  the  central  peak  due  to  structural  relaxa- 
tion has  a close  analogy  to  the  dynamic  central  peak  associated  with 

12 

the  structural  phase  transition  in  ferroelectric  crystals  as  well 
as  with  the  central  peak  in  orientationally  disordered  crystals 
which  arise  as  a result  of  the  dynamic  coupling  of  translation 
and  rotation.  In  fact  the  treatment  given  in  the  present  work  is 
closely  parallel  to  that  given  in  ref.  13  for  dynamics  of  translations 
and  rotation  in  molecular  crystals. 


C)  Rayleigh-Brillouin  spectral  density 

Using  Eq.  (29),  we  can  also  obtain  an  expression  for  the  Rayleigh- 
Brillouin  spectrum  as  a function  of  frequency.  For  this  we  separate 
A2(o>)  and  A 3 (ui)  into  real  and  imaginary  parts. 


A2(oi)  ■ a ^(oi)  + ia2"(u)) 


13(01)  * 013(01)  + iaj’Voi) 


(54) 

(55) 


Substituting  Eqs.  (54)  and  (55)  into  Eq.  (29)  and  after  carrying 
out  some  algebraic  manipulations,  we  obtain 
NjDi  + N2D2 


I (q.oi) 


Di2  + D22 


(56) 


where 


Ml  - - (1  + --— ^-)  (1  +— ^-)oi2  + ct' (oi)ai(oi)  + (1  -~)ui  2 (57) 


Y s 


«2  (w)  a 3 ' ' (m) 

&2  ' (ui)+a3  ' (oi)  + — at  3 ' (oi)  + a2'(o»)  V ' 


(58) 


* 


D2  = id 


_u2(1+  (1+^i)+a)2 


(1  + > + 0.2*  (ajjcta*  (a>) 


(60) 


The  detailed  forms  of  (^'(w),  oc2  ' ' (gj)  , ot  3 ’ (oj  ) and  ot  3 ’ ’ (oj  ) depend 

on  the  specific  time  dependence  of  the  memory  functions  K22(t)  and 

K33(t).  There  are  techniques  available  to  obtain  the  asymptotic 

frequency  dependence  of  A2(w)  and  A 3(0)).  For  example  if  short-time 

dynamic  processes  are  important  in  determining  the  memory  functions, 

one  may  expand  X ^ (to)  in  a power  series  in  w and  calculate  the 

14 

moments  of  memory  functions.  On  the  other  hand,  if  the  long-time 
dynamic  processes  are  important,  we  may  expand  A^Cu)  in  a power  series 
in  w.  However,  for  the  purpose  of  illustrating  the  effect  of  structural 
relaxation  it  suffices  to  use  here  the  functional  forms  given  by  Eq. 
(41)  for  the  calculation  of  the  Rayleigh-Brillouin  spectrum  in  the 
entire  frequency  region. 

In  the  next  section,  we  shall  apply  the  present  theoretical  re- 
sults to  interpret  the  Rayleigh-Brillouin  spectrum  of  a highly  viscous 
polymeric  liquid  such  as  polypropylene  glycol.  We  shall  show  that  the 
structural  relaxation  plays  a very  important  role  in  determining  I(q,w) 
and  in  order  to  account  for  the  experimental  results  it  is  very 
important  to  retain  the  frequency  dependence  in  A 2 (to) , although  the 
frequency  dependence  of  A3(u>)  is  unimportant.  In  the  subsequent 
section  we  supress  the  frequency  dependence  of  A3  and  equate  it  to 
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IV.  Comparison  With  Experimental  Results 

We  now  use  Eq.  (56)  together  with  Eqs.  ( 57 ) — ( 60)  and  Eq . (41) 
and  a least  squares  fitting  procedure  to  calculate  the  Rayleigh- 
Brillouin  spectrum,  using  the  parameters  suitable  for  polypropylene 
glycol  (PPG).  The  data  considered  are  the  frequency  shift,  linewidth 
and  the  Landau-Placzek  ratio  as  a function  of  .emperature. 

For  the  comparison  we  first  use  a trial  and  error  procedure  to 
compute  the  spectrum  at  high  temperature  such  that  the  theoretical 
spectrum  fits  the  measured  one.  This  gives  a set  of  reference  data 
points  needed  for  subsequent  calculations  of  the  theoretical  spectra 
at  other  temperatures. 

The  reference  data  set  consists  of  uig,  t2,  cij  and  appropriate 
for  the  high  temperature  (400  K)  spectrum  of  PPG  M.W.  425.  To 
obtain  results  for  the  spectra  at  other  temperatures,  we  assume  that 
t 2 has  an  Arrhenius  temperature  dependence 

t2  =*  t2°  exp(Ea/kT)  (61) 

This  assumption  is  justified  from  the  experimental  results  for 
4 

PPG.  Thus  instead  of  t2  we  replace  it  by  0 and  Ea  for  calculating 
the  results  at  other  temperatures. 

Since  accurate  experimental  data  for  the  frequency  shift  (a)g) , 
linewidth  (O  and  intensity  are  available  for  PPG-425,1  we  have  found 
that  the  following  least  squares  fitting  procedure  is  most  economi- 
cal in  computer  time.  We  used  Eq.  (56)  to  obtain  theoretical  values 
for  dig  and  T and  then  compared  with  the  observed  tug  and  T until  the 
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standard  deviation  (in  the  least  squares  sense)  between  the  experi- 
mental and  theoretical  values  are  simultaneously  minimized.  This 
procedure  gives  a tentative  set  of  parameters  w , T2°,  E , a 2 and 
D^,  which  are  consistent  with  the  observed  frequency  and  linewidth 
values.  To  improve  the  theoretical  fit,  we  calculate  the  intensity 
and  compare  the  results  with  the  experimental  data.  This  is  accom- 
lished  by  a fitting  cycle  to  check  the  self-consistency  of  all 
experimental  data  available.  In  cycle  1,  we  use  the  best  values 
parameters  u)g,  t2°,  Ea,  a 2 and  D^.  This  set  gives  a good  fit  to  all 
values  of  w and  T,  but  it  may  not  fit  the  intensity  properly.  The 

O 

set  of  parameters  is  then  varied  to  generate  the  best  integrated  inten- 
sity results  for  Rayleigh  and  Brillouin  peaks.  In  this  process,  a new 
set  of  parameters  are  generated  to  fit  the  frequency  and  linewidth. 

The  cycle  is  repeated  until  the  set  of  parameters  fits  simultaneously 
all  of  the  uig,  T and  intensity  results. 

In  the  process  of  curve  fitting,  we  have  found  that  Che  intensity 
of  the  central  Rayleigh  component  is  affected  significantly  by  the 


value  of  Dt-  Furthermore,  we  have  also  found  that  the  values  obtained 
for  can  be  correlated  with  an  Arrhenius  equation  given  by 

Dt  - Dt°  e“Et  /VtT  (62) 

with  D^,°  » 0.98(GHz)/q2  and  Et°  ■ 3.16  kcal/mole.  This  gives  the  value 
of  yD^fl2equal  to  4.9  MHz  at  300  K.  The  theoretical  fits  to  the  experi- 
mental data  of  to  and  T are  also  improved  if  the  temperature  dependence 

O 

of  wg  is  allowed.  This  is  expected  as  the  adiabatic  sound  frequency  is 
expected  to  decrease  with  increasing  temperature.  Thus  for  convenience 


we  assume 


m - u>  ° ( 1-mT) 

9 S 


(63) 
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The  best  values  of  the  parameters  are  j ° = 11  GHz  ana  m = 0.0017  K 

s 

o - l 4 

Using  these  we  obtain  E^  = 4.4  kcal/mole,  t2  = 2.36x10  sec  and 
a2  " 0.438  for  the  best  fit.  a,  is  fairly  insensitive  to  the 
temperature  change.  We  have  tried  other  temperature  dependences  for 
w and  have  also  obtained  good  fits  to  both  w_  and  r , but  the  results 

^ D 

obtained  for  E differ  within  15%  at  the  most.  This  is  about  the 

a 

experimental  error  involved  in  determining  E^.  With  all  of  the  para- 
meters available,  we  have  calculated  the  m and  P values  as  a function 
of  temperature.  The  results  are  shown  in  Fig.  1,  together  with  the 
experimental  values  for  comparison.  The  agreement  is  quite  good.  The 
slight  deviation  for  the  range  between  300  and  350  K is  due  to  the 
assumed  simple  form  of  Eq.  (63).  The  calculated  Landau-Placzek  values 
for  PPG-425  for  the  temperature  range  between  185  and  400  K are  shown 
in  Fig.  2,  together  with  the  experimental  results. 

Due  to  the  presence  of  the  structural  relaxation,  a few  words 
about  the  determination  of  the  Landau-Placzek  ratio  are  in  order.  In 
Fig.  2,  the  calculated  intensity  ratio  excludes  the  intensity  contri- 
bution from  the  structural  central  peak.  While  this  is  easily  done 
theoretically,  it  is  difficult  to  accomplish  experimentally,  especially 
in  the  temperature  region  where  the  structural  relaxation  is  most 
pronounced.  Efforts  were  made  to  minimize  the  contribution  from  the 
structural  relaxation  in  determining  the  accurate  Landau-Placzek  ratio. 
This  was  done  by  first  integrating  the  Brillouin  intensity  beyond 
the  peak  frequency.  The  total  Brillouin  intensity  of  each  peak  is 
assumed  to  be  twice  this  value.  This  thus  neglects  a small  correction 
due  to  lineshape  asymmetry,  and  due  to  the  intensity  at  the  wing  from 


i 


the  structural  central  peak.  The  error  arising  from  this  procedure 
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is  insignificant  at  both  high  and  low  temperatures.  However,  in  the 
temperature  range  where  the  structural  relaxation  is  most  significant 
the  error  is  somewhat  larger.  The  Rayleigh  intensity  was  determined 
by  first  extrapolating  the  wing  of  the  Rayleigh  peak  to  the  baseline 
and  then  integrating  the  intensity.  This  tends  to  over-estimate  the 
intensity  by  including  a small  contribution  due  to  the  structural 
central  peak,  but  the  error  is  believed  to  be  small  in  this  case.  The 
Landau-Placzek  ratio  data  determined  this  wav  are  thus  more  reliable 
than  those  reported  in  ref.  1.  One  notes  in  Fig.  2 that  the  agreement 
between  theoretical  and  experimental  results  is  rather  good.  The 
larger  discrepancy  in  the  temperature  range  between  220  and  300  K is 
presumably  due  to  the  experimental  difficulty  in  separating  accurately 
the  Rayleigh  and  the  Brillouin  intensities  from  the  complete  experi- 
mental spectra. 

Using  the  obtained  best  set  of  parameters  appropriate  for 
PPG-425,  we  have  then  computed  the  Brillouin-Rayleigh  spectra  using 
Eq.  (56),  together  with  Eqs.  (57)— (60 ) and  Eqs.  (49)  and  (50). 

Six  representative  theoretical  spectra  at  400,  350,  310,  280,  240 
and  200  K are  shown  in  Fig.  3A-F.  These  spectra  clearly  demonstrate 
the  evolution  of  the  effect  of  structural  relaxation  on  the  Rayleigh- 
Brillouin  spectrum.  These  spectra  when  convoluted  with  the  instru- 
mental function  are  in  good  agreement  with  the  observed  spectra,  thus 
indicating  that  the  present  theory  is  adequate  in  interpreting  the 
Rayleigh-Brillouin  scattering  of  polymeric  liquids  in  which  the 
structural  relaxation  plays  a very  important  role.  In  Fig.  3,  one 
notes  that  while  the  width  of  the  Brillouin  doublet  increases  with 
decreasing  temperature  between  400  and  310  K,  there  is  an  increase 
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in  the  intensity  in  the  region  between  the  Ravleigh  and  the  Brillouin 
peaks.  This  is  associated  with  the  central  peak  due  to  structural 
relaxation,  similar  to  the  Mountain  peak  observed  in  simple  molecular 
liquids.^’  ^ However, as  the  temperature  is  decreased  further,  the 
Brillouin  peaks  narrow  and  the  structural  central  peak  becomes 
gradually  localized  to  the  wing  area  of  the  thermal  diffusion  Rayleigh 
peak.  The  behavior  of  the  structural  central  peak  is  remarkably 
sensitive  to  the  temperature  variation. 

V.  Summary  and  Conclusion 

We  have  used  the  linear  response  theory  to  analyze  the  Rayleigh- 
Brillouin  scattering  spectrum  of  a viscoelastic  liquid.  The  theory 
includes  a complete  set  of  dynamic  variables  involving  density, 
velocity  and  energy  fluctuations.  An  exact  equation  including  the 
effect  of  structural  relaxation  has  been  obtained.  In  the  fast 
relaxation  limit  it  is  shown  that  the  Rayleigh-Brillouin  scattering 
is  identical  to  that  predicted  by  the  classical  hydrodynamic  equation. 
New  results  are  obtained  for  the  slow  relaxation  case.  In  this  case 
it  is  shown  that  in  addition  to  the  unshifted  Rayleigh  peak  a new 
central  peak  is  also  present . This  new  peak  is  due  to  dynamic 
coupling  of  the  translational  motions  of  the  scattering  particles 
with  collective  structural  relaxation,  in  analogy  to  the  Mountain 
peak  which  arises  from  the  coupling  of  translation  to  intramolecular 
vibration  of  simple  molecules.  Our  theoretical  results  for  the 
Brillouin  peak  and  the  central  peak  linewidth  differ  from  those  ob- 
tained by  Mountain  in  the  case  of  thermal  relaxation.  Mountain's 
calculation  is  based  upon  modified  hydrodynamic  equations  and  ours  is 
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in  terms  of  the  generalized  Langevin  equation  derived  by  using  Zwanzig- 
Mori's  projection  operator  technique.  The  structural  central  peak 
considered  here  is  analogous  to  the  dynamic  central  peak  associated 
with  the  phase  transition  in  solids. 

We  have  described  the  evolution  of  structural  relaxation  and  its 
effect  on  the  Rayleigh-Brillouin  spectrum.  The  structural  central 
peak  is  most  pronounced  when  the  frequency  dispersion  and  the  line- 
width  maximum  are  present.  We  have  also  used  the  theoretical  result 
to  calculate  the  Brillouin  frequency  shift,  spectral  linewidth  and 
the  Landau-Placzek  ratio  as  well  as  the  entire  Rayleigh  Brillouin 
spectrum  as  a function  of  temperature.  We  show  that  the  theory 
including  structural  relaxation  gives  a good  fit  to  the  experimental 
results  of  polypropylene  glycol. 
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Appendix 


We  show  here  Chat  the  dynamic  variable  £3  defined  in  Eq.  ( 8 ) 
in  the  text  is  proportional  to  the  temperature  fluctuation. 

Since  5p(q)  and  5e(q)  are  related  to  5p(r)  and  <5e(r)  by^ 


6p(q)  - [ 5p(r)ei.  - d3r 


(A—  1 ) 


5e(q) 


6e(r)ei3’5  d3r 


(A-2) 


where  5p(r)  and  5e(r)  are  respectively  the  fluctuations  in  number 
density  and  energy  density. 

If  the  system  under  consideration  is  in  local  equilibrium,  then 
from  the  first  law  of  thermodynamics,  we  have 


5e  «*  T <$ .s  + ( — — — ) 5p 
0 po 


(A- 3) 


where  T , p , e and  p are  respectively  the  temperature,  number 
0000 

density,  energy  density  and  pressure  at  equilibrium.  5s  * 5S/V, 

V being  the  volume  and  5s  the  entropy  fluctuation  per  unit  volume. 
Using  thermodynamics,  one  also  has 


mp  C 

a 0 v 

4,1  + "I—  5T 


(A-4) 


5p  and  5T  are  independent  thermodynamic  variables  so  that  <5p5T>  • 0. 
Substituting  Eq.  (A-4)  into  Eq.  (A-3),  we  obtain 


(eo+  M ii-1 

■'  °o  ' poxT. 


5p  + mp^C  6T 
0 v 


(A- 5) 
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and 


<6e5p>  • 


c + D I 01 T 


, t u > 

(^f-y 


°0X1 


< (5p) 2> 


(A-6) 


After  establishing  Eq.  (A-6),  we  consider  £3  (Eq.  (8)) 


( /l  - n2  ) <|6e(q)|2>ii 


5e(q)  - 


:|6o|2> 


(A-7) 


<5E(q)<Sp*(q)>6p(q) 

Using  Eqs.  (A— 1 ) and  (A-2) , we  can  write 

<6e(q)6p*(q)> 

d3rd3r'  e1?  ^5  5 ^ <<5e(r)6p  (r ' )> 

d3rdV  eiq,<^,) 


'o+po 

aT 

0 

f 

°0 

" v . 

J 

VPo 

oIo  1 

< 1 A 

. °0 

X 

0 

a. 

<i  i 

: | 5p (q) 1 : 


Substitution  of  Eq.  (A-8)  into  Eq.  (A-7)  gives 


<<5p(r)5p(r')> 


(A-8) 


53  - 


1 


< 1 6 e ( q ) | 2 >'*  /I  - n2 
6p(q) 


5e(q) 


/£o+Po  aTo  \ 
\ po  poXT  / 


(A-9) 

Using  Eq.  (A-5) , one  can  reduce  this  equation  immediately  to 
mp  C 

*3 3 ~ ' — 6T(q)  (A-10) 


< | <$e  | 2>^  /T-  n2 


Thus  £3  is  proportional  to  <5T(q). 
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Figure  1. 


igure  2. 


Figure  Captions 

Comparison  of  Che  measured  and  calculated  Brillouin 
peak  frequency  and  linewidth  data  of  PPG  425  as  a 
function  of  temperature. 

Comparison  of  the  measured  and  calculated  Landau- 
Placzek  ratio  data  of  PPG  425  as  a function  of  tem- 
perature . 

Simulated  Rayleigh-Brillouin  spectra  for  PPG  425. 

(A)  at  400  K,  t 2<i>s  - 0.12;  (B)  at  350  K,  t2ws  “ 0*34; 

(C)  at  310  K,  T2wg  - 0.91;  (D)  at  280  K,  i2u>g  - 2.2; 

(E)  at  240  K,  t2w  *9.1;  and  (F)  at  200  K,  t2oi  - 64. 
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